Multiaxial Kitagawa Analysis of A356-T6 



M. J. Roy a , Y. Nadot b >*, C. Nadot-Martin b , P.-G. Bardin b , D. M. Maijer a 

a Dept. of Materials Engineering, The University of British Columbia, Vancouver, BC, 

Canada V6T 1Z4 

b Institut PPRIME - CNRS - Universite de Poitiers - ENSMA - UPR 3346 - 
Departement Mecanique des Materiaux - Teleport 2-1 Avenue Clement Ader - BP 
40109 - 86961 FUTUROSCOPE CHASSENEUIL Cedex - France 



Abstract 

Experimental Kitagawa analysis has been performed on A356-T6 contain- 
ing natural and artificial defects. Results are obtained with a load ratio of 
R — — 1 for three different loadings: tension, torsion and combined tension- 
torsion. The critical defect size determined is 400 ±100 /im in A356-T6 
under multiaxial loading. Below this value, the microstructure governs the 
endurance limit mainly through Secondary Dendrite Arm Spacing (SDAS). 
Four theoretical approaches are used to simulate the endurance limit charac- 
terized by a Kitagawa relationship are compared: Murakami relationships [Y 
Murakami, Metal Fatigue: Effects of Small Defects and Nonmetallic Inclu- 
sions, Elsevier, 2002.], defect-crack equivalency via Linear Elastic Fracture 
Mechanics (LEFM), the Critical Distance Method (CDM) proposed by Sus- 
mel and Taylor [L. Susmel, D. Taylor. Eng. Fract. Mech. 75 (2008) 15.] and 
the gradient approach proposed by Nadot [Y. Nadot, T. Billaudeau. Eng. 
Fract. Mech. 73 (2006) 1.]. It is shown that the CDM and gradient methods 
are accurate; however fatigue data for three loading conditions is necessary 
to allow accurate identification of an endurance limit. 

Keywords: A356 T6, Eshelby inclusion, Casting defect, Multiaxial 
Kitagawa diagram, Critical defect size 
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1. Nomenclature 

v^area Defect size parameter defined as the square root of a defect cross- 
sectional area 

v /area ref Reference defect size 

b g Gradient model material parameter 

/_i Endurance limit under fully reversed tension 

a,A,k Murakami parameters 

t-i Fatigue limit under fully reversed torsion 

F{b/a) Function of defect aspect ratio (b/a) 

H v Vickers hardness 

Ii t a First invariant of the stress amplitude tensor 

Ii >m First invariant of the mean stress tensor 

J 2 Second invariant of S a 

S_ Deviatoric stress tensor 

</2,mean Mean value of J 2 over a period 

J 2 max Maximum value of J' 2 over a period 

R Load ratio where R = o- min / cr max 

Y Crack shape factor 

7i> 72, 73, P Material parameters 

E Stress tensor 

a max Maximum stress over a loading cycle 
o"min Minimum stress over a loading cycle 
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a c Critical stress amplitude 
<j\ Maximum principal stress 
0"2 Minimum principal stress 

aw Endurance limit of a specimen with a defect under tensile loading 

t w Endurance limit of a specimen with a defect under torsional loading 

cr eqvu , cr* q Equivalent stress 

cr eqvu Max Defect tip equivalent stress 

cr eqvu Equivalent far-field stress 

Acr c LEFM critical stress amplitude 

A-K"th,eff Effective stress intensity factor threshold 

Vcr eqvu Equivalent stress gradient 

T Period length 

t Time 

(Tf e Experimental endurance limit 
<7f c Calculated endurance limit 

5 C Percent difference between calculated and experimental endurance limits 
such that 5 C = (cr/ e — er/ c )/er/ e x 100% 

2. Introduction 

The tensile fatigue behaviour of A356-T6 has been the subject of study 
by a number of researchers [H, Q, S 0, 0, @, 0]. In almost all related stud- 
ies, casting defects such as intermetallic inclusions, porosity, shrinkages and 
oxide films have been shown to be present at the origin of the failure. In 
cases where defects were not present, the basic microstructure jil, Isf has been 
shown to determine the fatigue life. The presiding microstructural factor 
in this latter case is Secondary Dendrite Arm Spacing (SDAS) which de- 
cides the overall strength of the material. By processing samples of A356 
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via Hot Isostatic Pressing (HIP) to eliminate porosity, Gao et al. [8j com- 
pared the influence of SDAS and porosity on the tensile endurance limit. It 
was found that the fatigue limit was significantly increased either through 
HIP or halving the SDAS. While the effects of porosity and SDAS are not 
mutually exclusive, the role SDAS plays is less important when the material 
contains defects. Another microstructural characteristic that participates in 
fatigue mechanisms is the secondary eutectic phase. Fatigue crack initiation 
has been found to occur at silicon particles within this phase jE), in sam- 
ples that were free of defects. However, in the majority of prior studies, the 
authors did not explicitly quantify the critical defect size. One exception is 
the work of Brochu et al. |5| where an experimental Kitagawa relationship 
was developed for rheocast A357, a similar alloy to A356. This study demon- 
strated a critical defect size of 150 /im under fully reversed (R = —1) tensile 
loading. In terms of multiaxial fatigue behaviour, only two previous studies 
are available for A356 0, [To|. De-Feng et al. 10] performed tension-torsion 
fatigue testing on thin-wall specimens in the low-cycle regime (10 4 cycles). 
McDowell et al. jsf performed torsional High-Cycle Fatigue (HCF) testing, 
however these tests were conducted with deformation control. 

The fatigue life prediction of nominally defective materials such as A356 
is of great importance to industry and has been the subject of considerable 
study. Linear Elastic Fracture Mechanics (LEFM) has been shown to only 
apply to the study of long cracks where stress fields are homogeneous and not 
affected by local plasticity, making it inappropriate to be applied on short 
cracks ll|. Therefore, treating defects as cracks and leveraging analytical 
approaches based on LEFM may not always be appropriate. Using experi- 
mental assessments of local plasticity through microhardness measurements, 
Murakami 12] has shown that fatigue behavior can be assessed for some ma- 
terials through empirical relations based on applied stress and hardness. In 
terms of more advanced computational and analytical assessments, there are 



two stress-based approaches as highlighted by Atzori et al. [13[ . The first are 
local methods which dictate that an effective stress must be reached around 
a defect to affect the fatigue resilience 



14, 15 or nominal stress methods 



where the fatigue or endurance limit of the material is defined in terms of a 
nominal applied stress 16j. Most recently, Critical Distance Theory 17| has 



been shown to successfully correlate fatigue behavior with a critical distance 
from a defect that is material dependent. Other studies have shown that 
the stress gradient around a defect is a more practical indicator of fatigue 



resilience 18, 19, 20 
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The objective of the present study is to investigate the influence of cast- 
ing defects on the HCF behaviour of A356-T6. Kitagawa-type analysis is 
performed with experimental results for three different scenarios: tension, 
torsion and combined loading. The critical defect size from a multiaxial 
standpoint is then defined. Finally, four different approaches to simulating 
the evolution of the endurance limit with increasing defect sizes are compared 
for each of the loading cases. This is done to determine the best method for 
simulating the endurance limit based on defect size for this industrially rel- 
evant alloy. 

3. Material and experimental conditions 

The material employed in this study was Low-Pressure Die Cast (LPDC), 
strontium modified A356 (Al-7Si-0.3Mg) in the T6 condition with a typical 
chemical composition given in Table [TJ The majority of fatigue specimens 
were cut from a wedge-shaped casting, and a lesser number were cut directly 
from an automotive wheel casting. While both casting types were made 
with permanent steel dies, the wheel casting was actively cooled during so- 
lidification while the wedge casting was left to cool passively. As these two 
casting types encompass a wide range of solidification conditions, so too did 
the specimens from a defect and microstructural standpoint. The fatigue 
behaviour characterized in the current work is thus expected to represent 
material produced from a range of industrial practices. 



Tabic 1: A356 composition in wt-% 



Element Si Mg Na Sr 

Range (wt-%) 6.5-7.5 0.25-0.4 -0.002 -0.005 



3.1. Material preparation and micro structure 



A wedge was generated via the steel die, shown in Figure 1(a) The cast 



ing produced was 250 mm in height, 300 mm in length and had a 7° wedge 
angle. The motivation for this geometry was to create a gradient in cooling 
rate varying with height in the wedge. A more refined microstructure was 
expected at the base of the wedge where the cooling rate was the highest 
as compared to a coarser microstruture at the top where the cooling rate 
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was the lowest. One half of the wedge was devoted to tension-torsion fa- 
tigue specimens only, and the other half had a mixture of tension, torsion, 
and combined tension-torsion fatigue specimens (Figure 1(b) ) drawn from it. 
All specimens were then heat treated to a T6 condition after being removed 
from the wedge block with the following schedule: solutionized at 538°C for 3 
hours, quenched in water at 60°C, and finally artificially aged at 150°C for 3 
hours. Tension-torsion specimens were extracted from the spokes of the com- 
mercially cast automotive wheel which had the same material composition 
as the wedge casting. The wheel was subjected to the same T6 treatment 
before specimens were extracted. Tensile testing of wheel material resulted 
in a modulus of elasticity of 66 GPa, a yield strength of 164 MPa and an 
ultimate tensile strength of 317 MPa. 

Hemispherical artificial defects were introduced to the middle of the gage 
sections of six specimens via Electro-Discharge Machining (EDM). In total, 
four tensile and two torsion specimens were drawn from the bottom of the 
wedge and had artificial defects applied post heat-treatment. This technique 
of generating artificial defects has been qualified in other crack propagation 



investigations 13, |22|, and the size of these defects is presented in Table [3J 



A cross-section of a typical artificial defect is given in Figure 3(c^ 



Metallographic specimens were cut from the wedge and wheel at loca- 
tions corresponding to where fatigue specimens were taken. This resulted in 
samples from 8 evenly spaced locations through the height of the wedge, and 
3 locations corresponding to the wheel specimens. On each of these samples, 
the SDAS and porosity were characterized via optical microscopy and Clemex 
Vision PE software. Porosity was assessed by both the area percent and the 
y / area of each pore. Specimens and metallography results were grouped to 
provide four families of specimens: one for the wheel, and one each for the 
low, mid-range and high SDAS wedge specimens (Table [2]). The average 
macrohardness value H v value was found for each specimen family with a 
load of 5 kg and a minimum of 30 indents on two specimens per specimen 
family. The average H v value for all material was 85 MPa. Both the A/area 
parameter of all defects and the overall hardness of this material are within 
the range as specified by Murakami 12 



The secondary dendrite arm spacing showed a definite increase from the 
bottom of the wedge to the top coinciding with the cooling rate differential 
imposed by the casting practice. Average porosity remained relatively uni- 
form throughout the wedge based on the percent area, with the peak pore 
size found through metallographic analysis was 94 /xm. These porosity mea- 
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(a) Wedge mold and the resultant wedge 
casting 




(b) The three different types of fatigue 
specimens and their respective geome- 
tries: (i) tension-torsion, (ii) torsion and 
(hi) tension type. 

Figure 1: Wedge mold, resultant casting and the different types of fatigue specimens used 
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surements show that the middle of the wedge had higher porosity than the 
top and bottom of the wedge as characterized by the area percent and max- 
imum pore size. Furthermore, while the bottom of the wedge approached 
the same SDAS as the wheel, the porosity in the wedge was approximately 
double that of the wheel casting. 



SDAS Porosity H v ±1 St. Dev. 

Family (/^ m ) Area Max. ^/area (MPa) 

(%) (H 

Wheel (W) 37 ±8 0.0603 52 86.15 ±2 

Wedge Bottom (A/B) 42 ±11 0.1237 32 78.9 ±5 

Wedge Middle (M) 59 ±21 0.1244 94 83.1 ±3 

Wedge Top (T) 68 ±26 0.1232 42 89 ±7 

Table 2: Secondary dendrite arm spacing, porosity measurements and Vickers hardness 
values for all material. Family A denotes specimens extracted from the bottom of the 
wedge which had hemispherical artificial defects applied via EDM. Family B denotes spec- 
imens extracted from the bottom of the wedge containing only natural defects. 



3.2. Fatigue testing 

All fatigue tests were run in load controlled mode with a sinusoidal sig- 
nal and fully reversed (R = —1). The combined tension-torsion specimens 
were tested on an Instron servo- hydraulic test platform at 11 Hz, while the 
pure torsion and pure tension specimens were tested at 45 Hz on a Amsler- 
Vibraphore machine, atigue testing was conducted by employing the step 



technique as originally outlined by Maxwell and Nicholas [2l|. Each sped 



men underwent cyclic loading starting at a given stress amplitude assumed 
to be below the expected endurance limit, and samples that did not fail af- 
ter 10 6 cycles were cycled again at a stress amplitude one step higher. This 
was repeated until the specimen ruptured after less than 10 6 cycles. The 
endurance limit (07J was taken to be the applied stress amplitude (<r a , r a ) 
of the final step. 

For example, specimen Wl underwent 10 6 cycles at r a =80 MPa without 
failure. The stress amplitude was increased by 5 MPa to r a =85 MPa, where 
it also withstood 10 6 cycles. The stress amplitude was increased again by 
5 MPa to r a =90 MPa and the sample failed after 722,000 cycles, having 
withstood 3 loading steps of 5 MPa per step. Table 3 describes the details 
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Table 3: Test history for all fatigue specimens. Specimens Al through A6 contained 
hemispherical defects applied via EDM. 



Specimen 


Loading 


; at Failure (MPa) 




Step 


Nf 


\/area 


Name 


Type a 






Numb 


er MPa/step 


(xlO 5 ) 


(pm) 


Wl 


TT 





90 


3 


5 


7.22 


59 c 


W2 


TT 





85 


l b 


N/A 


3.00 


59 c 


W3 


TT 


70 


70 


2 


5 


1.05 


59 c 


Bl 


TT 


70 


70 


2 


5 


0.76 


90 c 


B2 


To 





70 


2 


5 


3.98 


39 c 


B3 


TT 





100 


3 


10 


1.51 


30 c 


B4 


TT 





110 


2 


10 


8.83 


38 c 


Ml 


TT 


95 





3 


5 


0.79 


90 c 


M2 


TT 


65 


65 


2 


5 


2.46 


514 


M3 


TT 


70 


70 


3 


5 


5.26 


53 c 


M4 


To 





60 


2 


10 


3.25 


531 


M5 


To 





55 


l b 


N/A 


2.27 


90 c 


Tl 


TT 


65 


65 


5 


5 


4.24 


112 c 


T2 


TT 


65 


65 


2 


5 


1.29 


265 


T3 


TT 


65 


65 


1 


N/A 


9.08 


300 


T4 


TT 


60 


60 


l b 


N/A 


4.05 


496 


T5 


To 





50 


2 


10 


4.84 


265 


T6 


Te 


90 





5 


10 


6.63 


372 


T7 


To 





50 


l b 


N/A 


7.33 


310 


Al 


Te 


90 





l b 


N/A 


4.24 


398 


A2 


Te 


90 





3 


10 


1.29 


514 


A3 


Te 


80 





4 


10 


9.08 


740 


A4 


Te 


70 





2 


10 


5.26 


760 


A5 


To 





70 


4 


10 


4.84 


465 


A6 


To 





50 


2 


10 


6.63 


708 



a TT: tension-torsion, To: torsion, Te: tension 
b Failure before 10 6 cycles 
c Estimated based on the maximum ^/area found with metallography 



of the step testing testing conducted. While the results obtained from this 
testing are not endurance limits from a statistical standpoint, step testing 
is the only technique that permits the evaluation of an 'endurance limit' 
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for a natural defect of unknown size. Furthermore, since the material is 
non-ferrous and there were specimens with similar microstructure requiring 
more steps than others to rupture at the same stress amplitude, the effect of 
coaxing is considered negligible for this material. In the current work, the 
term 'endurance limit' is defined as the stress level at fracture for one million 
cycles. 

4. Experimental results 

Experimental fatigue test results are given in Figures El [4] and |6] in the 
form of bi-linear Kitagawa diagrams for each of the loading cases. The frac- 
ture surfaces were first examined macroscopically, and these observations 
guided Scanning Electron Microscopy (SEM) analysis of the initiation area. 
The goal of SEM observations (Figures El El and [7]) was to clearly identify the 
initiation site and where possible, measure the initiating defect. Multiaxial 
fracture surfaces were found to be tortuous and indentification of the defect 
at the origin of the fatal crack was non-trivial. The systematic methodology 
employed to find the defect and correctly identify the area of the defect on 
the Kitagawa diagrams presented is as follows: 

• Optical observation of the fracture surface employed to determine both 
the fatigue and final failure zones. 

• The fatigue zone of the fracture surface was found to be clear enough 
to identify a given initiation area. The 'river marks' observed on the 
surface converged to this unique zone. 

• If there was a clearly defined defect, then the size was assessed using 
the y / area parameter directly from the fracture surface. This was done 
regardless of the position of the defect relative to the free surface. 

For specimens where the initiation site was unidentifiable, the initiating fea- 
ture was estimated as the maximum -y/area of porosity found via metallog- 
raphy (Section 13. ip . The use of these results in the Kitigawa analysis is 
thus speculative. For specimens where there were multiple initiation sites, 
the largest identifiable defect closest to the surface was characterized. The 
results of this analysis for all specimens are summarized in the final column 
of Table El 
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Figure 2: Kitagawa diagram for A356-T6 under tension. Specimens Al through A4 had 
artificial defects applied. 



4-1. Tensile results 

Figure [2] presents the experimental Kitagawa relationship under pure ten- 
sion. In all samples, the initial defect size was readily identifiable on the 
fracture plane that was found to be perpendicular to the direction of the 
maximum principal stress. The primary finding from the tensile Kitagawa 
curve is that the critical defect size is relatively large: specimens T6, Al and 
A2 have a very low impact on the endurance limit (8 % reduction). These 
tensile specimens show that the material appears to be sensitive to defects 
only when a/ area is greater than than 500 /im. Figure [3] presents the fracture 
surfaces for specimens Ml, T6 and Al. Specimen Ml shows a oxide-related 
defect at the origin of the crack (Figure 3(a)) linking to subsurface porosity. 
Specimen T6 failed due to a 400 /im shrinkage pore that intersected with 
the surface of the sample (Figure 3(b)), while specimen Al failed due to a 
400 fim artificial defect (Figure 3(c)). Since specimens T6 and Al demon- 
strated the same endurance limit, it is concluded that the area parameter 
is able to correlate different types of defects, independent of the nature of 
the defect. Furthermore, as these two specimens exhibit the same endurance 
limit despite being at the upper and lower end of the range of SDAS, this 
defect size is used as the reference defect size (y / area ref ). Nevertheless, this 
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(a) 




x60 0BB0 EBftV 500JJFH 



(c) 

Figure 3: SEM observations of initiation sites on tensile specimens (a) Ml and (b) T6 
having natural defects; (c) Al having an artificial defect 

finding should be verified with larger defects having a greater impact on the 
endurance limit. In terms of artificial defects, fracture surfaces for specimens 
A2, A3 and A4 were very similar to Al, showing that the artificial defect 
was inarguably the initiation point. 

4-2. Torsion results 

The Kitagawa diagram for torsion testing, shown in Figure HI is similar to 
the tensile results (Figure |2]). The experimental points presented on the curve 
below 100 /im are for specimens that were unsuccessfully classified by frac- 
tography. The samples have been separated along the horizontal axis based 
on the porosity assessment of their location in the wedge to render individual 
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Figure 4: Kitagawa diagram for A356-T6 under torsion. Specimens A5 and A6 had 
artificial defects applied. 



tests identifiable from a stress amplitude standpoint. As the endurance limits 
vary from 55 to 95 MPa for specimens with unidentifiable defects, the Kita- 
gawa diagram under torsion exhibits a large amount of scatter as compared 
to the tensile results. The main complication in identifying defects in the 
sub 100 /im range is the tortuous nature of the fracture surface, as shown in 



Figure 5(a) Cracking was activated on two planes of maximum shear such 
that the final fracture surface reveals multiple initiation sites. While multiple 
initiation sites could explain the scatter seen in the Kitagawa plot, careful 
examination of the fracture surface at each suspected initiation point (Figure 



5(b) ) did not always result in an identifiable defect. An attempt to link the 



presence of porosity to the multiple initiation sites was made with specimens 
B2, M5, T5 and T7 by metallography performed on sectioned fracture sur- 
faces. There was little to no deviation found from the porosity measurements 
given in Table |2j In light of these findings, the critical defect size is difficult 
to assess. Specimens with identifiable defects show a definite decrease in 
endurance limit beyond 300 /im, which is smaller than under tension. As is 
the case with the tensile testing, this critical defect size should be clarified 
by other tests on samples containing 300 /im or larger artificial defects. 
Figure [5] presents an example fracture surface for a specimen under pure 
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(e) 

Figure 5: Observations of initiation site defects on torsion specimens (a,b) W2 (c,d) T7 
and (f) M4 
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torsion and gives an indication of the typical challenge experienced in de- 
termining the defect size. It is surmised that the crack lips contact during 
crack propagation, and an excess of oxide fouls the fracture surface masking 
both the defect size and type. Macroscopically speaking, this black oxide 
appears throughout the gauge section of the pure torsion specimens. Fig- 
ure 



is an 



5(a) shows the fracture surface of specimen W2 and Figure 5(b) 
image of a crack origination point on the same specimen which is absent of 
any identifiable defect. Specimen T7 showed multiple initiation points, but 



no defects were observed at the crack origins (Figure 5(c)). The defect at 



tributed to specimen T7 (Figure 5(d)) was an identifiable defect located on 



the fracture plane of the main crack. Specimen M4 provides an example of 



a large, readily identifiable defect at the origin of failure (Figure 5(e) ). 



4-3. Combined tension-torsion results 

The combined tension-torsion Kitagawa diagram, presented in Figure El 
includes results for specimens with natural defects only. The macroscopic 
fracture surfaces were similar to those of the tensile specimens: a flat surface 
in the plane perpendicular to the direction of the maximum principal stress 
with clear, identifiable initiation sites the exception of specimen Bl (Figure 



7(a)). For the specimens tested, the Kitagawa diagram shows a very small 



influence of a 500 /mi defect such as specimen M2 (Figure 7(b) ) exhibiting a 



large 500 /iin subsurface pore. Below this size, there is no apparent influence 



of defects on the endurance limit. Figure 7(c) reveals a shrinkage void at the 
origin of the failure on sample T3. It is of interest to highlight that specimen 
M2 has approximately the same endurance limit ( 67 MPa) as specimen T3, 
reinforcing the independence of defect type and dependence of defect size 
characterized by a/ area on the overall fatigue life. 

4-4- Experimental summary 

The experimental results demonstrate an average endurance limit of 90 
MPa for pure tension (/_i) and 80 MPa for pure torsion (t_i). Furthermore, 
these results clearly indicate that the defect size necessary to affect the en- 
durance limit is relatively large. The critical defect size ranges from 300 /im 
for pure torsion to 500 yum for pure tension and combined tension-torsion. 
This suggests such that an overall critical defect size necessary to diminish 
the multiaxial endurance limit is y/a.iea, re{ = 400 /im. Together with this 
reference defect size, the tension and torsion endurance limits will be used in 
identifying model parameters in the following sections. 
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Figure 6: Kitagawa diagram for A356-T6 under tension-torsion 



Another important qualitative observation is that when the defect is char- 
acterized by the -y/area parameter, an artificial defect diminishes the en- 
durance limit in the same manner as naturally occurring defects. Therefore, 
the i/area is a powerful method to describe the size influence of a defect. 
However, analyzing the complicated topology of the fracture surfaces cou- 
pled with multiple initiation sites precludes an accurate assessment of the 
size of critical defects. 

Within a Kitigawa framework, the pure torsion results show a great deal 
of scatter, and indicate the need for further investigation. In light of the 
difficulty in assessing critical defect size for some specimens, more experi- 
mentation is necessary to categorically determine the effects of small defects 
on the endurance limit of this material. In terms of large defects, as the 
combined tension-torsion results do not clearly show a drop in endurance 
limit as compared to the tension results, further assessment of this loading 
scenario and material is also desirable. 

5. Simulation of multiaxial Kitagawa relationships 

Four standard models to predict the endurance limit of defective materials 
were employed to simulate the behaviour of A356-T6. What follows is a 
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presentation of each model, followed by a comparison of the results of each 
model to experimental data in Section El The following models are able to 
predict the Kitagawa relationships for tension, torsion and combined loading: 

1. Linear Elastic Fracture Mechanics (LEFM) approach whereby a defect is 
considered equivalent to a crack; 



2. Murakami relationships [12| employing an empirical function of the Vick- 
ers hardness, load ratio and defect size; 

3. Critical Distance Method (CDM) based on the stress field around a defect 
proposed by Susmel and Taylor [13]; and 

4. Gradient approach based on the stress gradient surrounding the defect 



proposed by Nadot [19|, [20 



In the case of the CDM and gradient approaches, a multiaxial fatigue criterion 
and a description of the stress distribution around a defect is necessary. 
For the multiaxial fatigue criterion, the equivalent stress criterion recently 



proposed by Vu 23( will be used. This criterion determines the multiaxial 
behaviour for complex loadings using an invariant approach. The equivalent 
stress criterion as proposed by Vu is generalized to account for phase-shifted 
loading effects for both hard and soft metallic materials. This equivalent 
stress is given as: 



a cq Vu = \JllJ2if) 2 + 72^!,mean + l^f{h,a, h,m) < P (1) 

where J' 2 {t) invokes the stress amplitude tensor deviator: 



m = ]llm ■■ m = + ( 2 ) 

For fully reversed loading, J' 2 (f) = ^2 max - While J2,mean is evaluated over a 
full period T: 



T 



J: 



2, mean 



\ j 4(t)dt (3) 



o 



Without mean stress, J2,mean = 0. The first invariant function, If(I^ a , ii im ), 
is a linear function of the first invariants of the stress amplitude and mean 
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stress where 



If(h,a, h,m) = h,a + «A,m (4) 

Ii >a = -{maxii(t) - min/i(t)} (5) 

h, m = -{nuucii(t) + mm I^t)} (6) 

For fully reversed loading conditions, Ix, m = and // = Ii >a = T, xx +T lyy +T lzz . 
The values of j3 and 71-3 are material properties based on the strength of the 
material and were empirically identified by Vu with the torsional endurance 
limit t-i, the tensile endurance limit /_i and the ultimate tensile strength. 
As the ultimate tensile strength of A356-T6 is less than 750 MPa, 71 = 0.65. 
The torsional endurance limit at 10 6 cycles for A356-T6 is taken to be 80 
MPa, therefore t_i = (3 = 80 MPa. The final material coefficient, 73 is given 
as: 



73 = ; 6 — (7) 

J-l 

The second input data required by both the CDM and gradient ap- 
proaches is the stress distribution around the defect. To aid in this calcu- 
lation, the defect is considered a spherical void in an infinite, homogeneous, 
isotropic, elastic matrix subjected to uniform stress at infinity. Local stresses 



around the void are calculated using the equivalent inclusion method [24[ and 



solved analytically (25j. This analytical description of the elastic stress dis- 
tribution around a void shows excellent correlation with that described by 
the Finite Element Method. 



5. 1 . LEFM approach 

LEFM describes a crack propagation threshold through the amplitude 
of the stress intensity factor, which is a function of crack length and stress 
amplitude Aa c . The defect size characterized by the area parameter is trans- 
formed to an equivalent semi-circular crack. The relationship between the 
endurance limit and defect size is therefore given by: 
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where Y is the crack shape factor and AK t ^ e g is the effective stress intensity 
factor threshold for crack propagation. The values used for A356-T6 were 
Y = 2/tt and AA' t h, c fr = 1.5 MPay^rn. These were not experimentally de- 
termined, but were synthesized from a large compilation of published data 
[E 0, 27, 28, 29]. While the LEFM approach requires an experimentally 



determined effective threshold stress intensity factor and a description of the 
defect shape to estimate the crack shape factor, this approach can be used 
to account for multiaxial complex loading and account for load ratio effects. 

5.2. Murakami relationships 



Murakami [12J proposed the use of the a/ area parameter to describe 
the size of a surface defect. The basis for this representation were non- 
propagating crack observations within a small stressed region surrounding a 
defect. From this study, the endurance threshold supposedly corresponds to 
the crack growth threshold. It was also demonstrated that the maximum 
stress intensity factor if/ max is linearly related to the -y/area parameter for 
different crack geometry. Furthermore, it was found that the fatigue crack 
growth threshold could also be correlated to a given value of Vickers hard- 
ness. The result is an empirical equation based on the defect size (i/area) 
and macrohardness to predict the endurance limit of materials containing 
small defects: 

A(H, + 120) fl-RY 

where a = 0.226 + H v x 10~ 4 and A = 1.43 or 1.56 for surface and internal 
defects, respectively. For torsional loading and surface defects, the endurance 
limit is expressed as: 

0.93(^ + 120) fl-R\ a 

T W = \ la " (10) 

F(b/a) (v^r^) 1/6 V 2 J 

where F(b/a) = 0.8397 for spherical defects. For combined tension-torsion 
loading, the endurance limit is given by the following: 

_ A(H V + 120) fl-R\" 

ai+kai= (v^r \~) { ) 

with k = —0.18 for cracks emanating from a round defect. There are two pa- 
rameters that are required for the preceding relationships: the macroscopic 
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Vickers hardness (85 MPa) as well as the value of F(b/a), which is a constant 
for spherical defects. While requiring minimal parameters to be identified, 
the work of Murakami is limited due to the inability of the stress state de- 
scription to allow for general multiaxial loading. 

5.3. Critical Distance Method 



Based on the methodology proposed by Susmel and Taylor [17J, the main 
tenet of this approach is to describe the influence of the defect through com- 
putation of the equivalent fatigue stress at a given distance from a defect. 
This requires an evaluation of the stress field surrounding a defect as de- 
fined previously (Eq. [1]). As the equivalent stress is calculated at the critical 
distance L/2 from the tip of the defect, the criterion is expressed as: 

m axfa eqvu (^n =(3 (12) 



<r=a c 



2 



where a is the nominal applied stress. The critical distance L/2 is therefore 
the point where o" e q Vu = P- In order to identify the critical distance, the 
experimental endurance limit is employed for the case where Uooref = 85 
MPa, and A/area ref = 400 jttm (Section |4.4p . This results in a value of 79 /im 
for L/2. 

5.4- Gradient criterion 



Based on the initial work of Nadot 19J], Gadouini et al |20j described the 



impact of a defect on fatigue resilience by amplifying the equivalent stress by 
a function of the surrounding stress gradient. This criterion is expressed as: 

< = ^ Vu , Max ( 1 - &,-^M < P (13) 
with the equivalent stress gradient expressed as: 

Max ^ 

A/area 



^IVu.Max "eClVu.oo 

V(7 eq Vu = rZZZZ ( 14 ) 



where d e q Vu Max is the maximum equivalent stress calculated at the tip of 
the defect and cr oqVu is the nominal equivalent far-field stress. Having the 
same dimension as yj area, the length parameter b g requires evaluation before 
this criterion can be applied. This parameter was identified with the same 
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experimental data used for the CDM approach (Section I5.3p . with croo rc f = 85 
MPa and v /area ref = 400 /im such that: 



6. Comparison between simulations and experimental results 

Figures [8] M and [10] present the predictions made by each of the four 
approaches compared with experimental results for the three loading cases. 
Under tension (jSJ), the CDM and Gradient approaches describe the exper- 
imental results quite well. However, this is likely due to the single experi- 
mental point employed for identification of the necessary parameters. The 
trend of the fatigue limit versus defect size is also well described. Murakami's 
equation leads to non- conservative results, but the trend is well described. 
LEFM gives conservative results, with an acceptable endurance limit trend. 
Under torsion (jHJ), all approaches adequately describe both the fatigue limit 
and trend; however, Murakami relationship is again non-conservative and 
the CDM approach leads to a small allowable defect size. Under combined 
loading (fTDl . all approaches are non-conservative except Murakami. Again, 
the critical defect size is largely underestimated by the CDM. 

Table H] demonstrates the error between experimental and simulated en- 
durance limit results for selected conditions. The difference between the cal- 
culated and experimentally found endurance limits (5 C ) is presented in Figure 
[TTl Non-conservative assessments are identified by negative values and posi- 
tive values indicate conservative ones. The result obtained by averaging the 
absolute error given by each approach leads to the following result: LEFM = 
19 %, Murakami = 20 %, CDM = 11 % and Gradient = 9 %. Based on this 
comparative assessment, describing the effect of a defect through the elastic 
stress field (Gradient or CDM approaches) provides better results than the 
LEFM or Murakami methods. However, Murakami's equation remains very 
good at estimating the endurance limit of A356 for different defect sizes and 
loading cases using only the macroscopic hardness of the material. LEFM 
provides conservative estimates with approximately the same average error 
for the same identification cost: the only experimental parameter required 
is the effective threshold stress intensity factor for long cracks. Both the 
CDM and Gradient approaches are more accurate, however they require the 
determination of the elastic stress field around the inclusion as well as three 




(15) 
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Figure 8: Comparison between tensile simulations and experimental results. Circled points 
indicate artificial defects. 
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Figure 10: Comparison between tension-torsion simulations and experimental results 



experimental endurance limits, including one with a defect. It is important 
to note that computation of elastic stresses is quite important in these latter 
approaches for A356-T6 as the endurance limit is half the yield stress. This 
indicates that there is a very small amount of plasticity at the tip of the 
defect. 

7. Conclusion 

• For A356-T6 submitted to multiaxial fatigue loading, fatigue cracks 
can initiate either on casting defects or inside the microstructure. Both 
scales are in competition for the localization of cyclic plastic deforma- 
tion that leads to the initiation of a crack. 

• When a crack initiates at a defect, the various types of defect can be 
characterized by \J area: natural defects were found to have the same 
endurance limit as artificial defects. 

• The critical defect size has been found to be 400 ±100 /xm in A356-T6. 
This result was obtained for both artificial and natural defects under 
three fully reversed (R = —1) loading scenarios: tension, torsion and 
combined loading. 
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Tabic 4: Comparison between experimental and calculated endurance limits. Specimens 
Al, A2, A4 and A5 contained artificial defects. 



Specimen 


Loading 


LEFM 


5c (%) 
Murakami CDM 


Gradient 


T6 




18.8 


-21.9 


-2.35 


-3.53 


Al a 
A2 


Tension 


22.4 
31.8 


-18.8 
-14.1 



5.88 




8.24 


A4 




31.4 


-30.0 


-5.71 


-1.43 


A5 




6.15 


-16.9 


16.9 


15.4 


M4 


Torsion 


-3.64 


-36.4 





3.64 


A6 




-11.1 


-57.8 


-13.3 


-11.1 


T2 




9.23 


3.08 


24.6 


9.23 


T3 


Combined 


26.6 


4.69 


25.0 


15.6 


T4 




34.5 


-1.81 


21.8 


21.8 



a V / area ref condition 



• Further experimental effort is needed to better characterize the Kiti- 
gawa relationship in pure torsion, understand the effect of small defects 
on all loading scenarios and to study the impact of larger defects in 
combined tension-torsion loading. 

• Multiaxial Kitagawa-type relationships were simulated using four dif- 
ferent approaches: Murakami, LEFM, CDM and the Gradient method. 
Results show that Murakami relationships give mainly non-conservative 
results with an average error of 20%. LEFM provides mainly conser- 
vative estimates with an average error of 19%. The error in using the 
CDM and Gradient approaches are both equivalent with mainly con- 
servative results and an average error of 11 and 9%, respectively. 
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Figure 11: Quantitative comparison of experimental results for each of the four simulation 
approaches. Specimens Al, A2, A4 and A5 contained artificial defects. 
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